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I. INTRODUCTION 



ibrium in an optical cavity reveals a BEC even for 
^ in solid helium 4 He has stimulated the search of a 



The existence of a Bose-Einstein condensate (BEC) in an ideal 3D quantum Bose gas 
served as a useful physical concept in the study of theoretical models for superfluidity, 
which are usually associated with the presence of a BEC. It is interesting that the physical 
realization of a BEC in a dilute gas has given impetus to the study of BEC for its own 
sake with a view of understanding many-body systems In addition, the confinement of 
photons and molecules in thermal equi 
photons [2J. The existence of superflow 
BEC in solid helium thus establishing the existence of BECs in all three states of matter-gas, 
liquid, and solid. The emergence of a self-organized supersolid phase, both a superfluid with 
crystalline order simultaneously, formed by a BEC coupled to an optical cavity has been 
observed [4|. The quantum phase transition describing the supersolid, which is associated 
with a spontaneous broken spatial symmetry, is in quantitative agreement with the Dicke 
model of superradiance and is driven by an infinitely long-range interaction between the 
condensed atoms 0]. A nonlocal potential is found to favor a crystalline BEC for the ground 
state of two-dimensional interacting bosons Numerical techniques have been used to 
predict a novel supersolid phase for an ensemble of Rydberg atoms in the dipole-blocked 
regime confined to two dimensions, interacting via a repulsive dipole potential softened at 



short distances 



It is claimed that the superfluid droplet-crystal phase does not crucially 
depend on the dipolar form of the interaction at long distances Q]. 

In this work, we use the Bogoliubov inequality to establish the properties of the interpar- 
ticle potentials that allow for a BEC in crystalline order that can be regarded as a supersolid. 
It is shown that a local, two-particle potentials cannot result in a BEC with crystalline order 
in 2D and what is needed are nonlocal, two-body potentials of infinite range. 

II. INTERACTING BOSE GAS 

Consider the Hamiltonian for an interacting Bose gas 

H = J iri\r)(^-V 2 )i>(r) + J <ir**(r)V^(r)*(r) 

+ /*,/ *> / *, / irMnWMVln, r 2 , r 3 , (1) 
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where V ext (r) is a confining, external potential, V(ri, r 2 , r 3 , r 4 ) is the two-particle interaction 
potential, and ip(r) and $(r) are bosonic field operators that destroy or create a particle at 
spatial position r, respectively. The two-particle interaction potential V(ri, r 2 , r 3 , r 4 ) must 
satisfy the following general conditions: (i) translational invariance, (ii) Galilean invariance, 
(iii) identical particles, (iv) rotational invariance, (v) space-reflection invariance, (vi) time- 



reversal invariance, and (vii) hermiticity [7]. Therefore, in general, 

y(ri,r 2 ,r 3 ,r 4 ) = 5(r 1 + r 2 - r 3 - r 4 )(i*i - r 2 |V|r 3 - r 4 ), (2) 

which is referred to as a nonlocal potential. A mathematically simpler potential can be 
deduced from Eq. (2) if, in addition, 

(ri - r 2 |V|r 3 - r 4 ) = 5(r 1 - r 2 - r 3 + r 4 )V(|ri - r 2 |), (3) 

and so 

V(r 1 , r 2 , r 3 , r 4 ) = -Sfa - r 3 )5(r 2 - r 4 )V(|ri - r 2 |), (4) 

which is referred to clS db local potential. 

Macroscopic occupation in the single-particle state ip{v) result in the non- vanishing js] of 
the quasi-average ip(r) =< x£>(r) > and so the boson field operator 

$(r) = V>(r) + £(r), (5) 

with 



and 



El£ k '| 2 = l, (7) 
k' 

where iVo is the number of atoms in the condensate and V(D) is the D-dimensional "volume" 
and < 0(r) >= 0. The operator 0(r) has no Fourier components with momenta {k'} that are 
macroscopically occupied and so / drip'(r)i/)(r) = 0. The separation of ^(r) into two parts 
gives rise to the following (gauge invariance) symmetry breaking term in the Hamiltonian 

(i) 

H, ymm =/ rf r^( ri ) / d r 2 / rf r3/ rf r^(r 2 )[^(r 1 ,r 2 ,r3,r 4 ) + y(r 2 ,r 1 ,r3,r 1 ) W (r3)^(r 4 ) +ft .c. 
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= I dT 1 tf(r 1 )x(ri) + h.c. (8) 

The presence of this nonzero H sym m in the Hamiltonian gives rise to further macroscopic 
occupation in states other than the original state given by ip(r) and so the condensate 
wavefunction ^(r) gets modified by augmenting the single-particles states where macro- 
scopic occupation occurs. In such a case, macroscopic occupation in the state b would give 
rise to macroscopic occupation in the states a, such that a ^ b, whenever the matrix ele- 
ment < ab\V\bb > of the potential V, which is the last term in Eq. (1), does not vanish. 
For instance, macroscopic occupation only in the single-particle state with momentum p, 
which corresponds to a uniform condensate, does not give rise to macroscopic occupation 
in any other momentum state since the matrix element in the momentum representation 
< qp|V^|pp > vanishes by momentum conservation unless q = p. This consistency proviso 
requires that the correct condensate wavefunction if)(r) corresponds to that which gives rise 
to no symmetry breaking term in the Hamiltonian. That is to say, H symm vanishes for the 
correct condensate wavefunction ijj(r). 

For instance, macroscopic occupation in the single-particle states with momenta k, k ± 
q 1; k±q 2 gives rise, with the aid of the symmetry breaking term H symm and owing to linear 
momentum conservation, to additional macroscopic occupation in single-particle momenta 
states. Therefore, for ${r) to be orthogonal to both ip(r) and x(r) and so H syrnrn = 0, 
one must have macroscopic occupation in all the momentum states k + n 1 q 1 + n2q 2 , with 
ni, n 2 = 0, ±1, ±2, • • • . Accordingly, 

*< W = S £ «k + „, q , + ,« q2 e*+».<l.+^> ' = Ar„ k(r) . (9) 

V V / n\,ni=— oo 

Note that the BEC (9) generates a real-space, crystalline distribution of atoms since u-^(r) = 
u^(r + t m ) for any primitive lattice translation vector t m with e^'^ m = 1 for i — 1,2 and so 
(9) is of the Bloch form. For D = 2, the vector t m = m 1 a. + m2h, where m\ and m 2 can take 
all integer values and a and b are the edges of the unit cell, which form parallelograms given 
by the five Bravais lattices. Note that a • q 1 = 2tt, a • q 2 = 0, b ■ q 1 = 0, and b ■ q 2 = 2n. 
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III. BOGOLIUBOV INEQUALITY 



The absence or presence of a BEC in spatial dimensions D < 2 is based on Bogoliubov's 
inequality 

\({A,A}) > k B T\([C,A])\*/([[C,H],&]), (10) 

where H is the Hamiltonian (1) of the system, the brackets denote thermal averages, and 
the operators A and C are arbitrary provided all averages exist. Consider first the case that 
H has arbitrary local interparticle and external potentials, that is, V(ri, r2, r^, r^) is given 
by Eq. (4). 

Consider the following operators, 

C = J dre ik ' r ^ t (r)Vi(r) (11) 

and 

A = JdrJ dr'e- ik - r f{r)f*{r')ft{r)^(r'), (12) 
where k is arbitrary. Now, 

([AU]) = V(D)([(7,i]> = ^ (13) 

q 

([[C,A\,&]) = ^-N, (14) 

lib 

with 

^k = E^ +k = ^/dr|^(r)|V kr (15) 
k' 1 

with the aid of Eq. (6), where N is the total number of particles. Now, |A k | < 1 by the 
Cauchy-Schwarz inequality where the equality holds for k = 0, that is, Aq — 1 with the 
aid of Eq. (7). The vector q ^ {k'}, where {k'} is the set of condensate vectors for which 
f k , jL 0. For a BEC at rest, |£ k | 2 = |£_ k | 2 . Note that ^ q+k ^ for (q + k) G {k'} and 
^q+k = for q ^ {k'} and k G {k'}. 

It is important to remark that Eq. (14) holds for local interparticle potentials but need 
not hold true for the general case of nonlocal potentials [3] • Therefore, for a BEC to exist in 
spatial dimensions D < 2, the l/fc 2 -singularity that results from the Bogoliubov commutator 
in the denominator of Eq. (10) must me removed. For local interparticle potentials this can 
be accomplished only provided the origin k = is a limit point (or point of accumulation) 
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of condensates. This would correspond to both momenta q 1 and q 2 in Eq. (9) approaching 
zero. Of course, in this case the BEC (9) ceases to be of the Bloch form, albeit, remaining 
nonuniform. The case where the interparticle potential is nonlocal is discussed in Section 
IV. 

We sum the Bogoliubov inequality (10) over the single-particle momentum states in the 
set {k'} constituting the condensate, which includes an arbitrary neighborhood of the point 
of accumulation of the condensate at k' = that corresponds to a condensate at rest. We 
want to find an upper bound of the anticommutator {{A, A*}) = 2(AA jf ) + {[A\ A}). We 
extend the sum over the first term (AA^) over all values of k thus obtaining a larger upper 
bound 

MN NV 2 {D) > ^(AA f >, (16) 
k 

where the condensate wavefunction ipij) is orthogonal to the operator (p'(r) and we assume 
that the condensate density is bounded from above by |/(r)| 2 < M with 1 < M < oo 
since / <ir|/(r)| 2 = V(D). In (16) use has been made of the completeness relation for the 
momentum eigenstates and a negative term resulting from a single commutation has been 
dropped. Note that we are considering a condensate where all the single-particle states 
with momentum k' are occupied macroscopically with k' = a point of accumulation. In 
addition, we are supposing that the number of particles in the "volume" V(D) is fixed, that 
is, we are employing a canonical ensemble and so / drfo (r)tp(r) = ^k^k = ^ is actually 
the c- number N. 

Consider next the sum over k' of the commutator ([A*, A]), 

'E([A\A]) = N V 2 (D)'£ > [l-\A ls! \ t \, (17) 
k' k' 

with the aid of (13) and where £q + k = for q ^ {k'} and k G {k'}. This sum over the 
commutator is bounded from above provided the sum is restricted to values of k' that have 
a finite, upper bound. Now the right-hand side (RHS) of inequality (10) becomes 

* S T|<[*Mi>iv<tt*, = '^fYlum^L^ly (18) 

k' 

with the aid of Eqs. (13) and (14). Note that the RHS is bounded in the upper limit of the 
sum; however, it is the lower limit as k' — > for D < 2 where the sum may diverge which 
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would result in no BECs, viz., Nq = for T > 0. Combining Eqs. (16)— (18), we have for 
the Bogoliubov inequality, 

k k 

Therefore, the existence of a BEC for T > requires the convergence of the sum on the 
RHS of (19) over the macroscopically occupied single-particle momentum states k' of the 
condensate. Note that the sums in (19) over the condensate momenta can be approximated 
by integrals according to J2\r' J dk' and so 

The integral on the RHS has no infrared divergence since by (15), (1 — l^j^'l 2 ) vanishes 
quadratically as k' — > and so the l/fc 2 -singularity is removed thus allowing the existence 
of a BEC for D< 2. 



IV. SUPERSOLID BEC 



The analysis of the possible existence of BEC for D < 2 in the previous section was based 
on two-particle, local interactions. Clearly, local interparticle potentials cannot give rise to 
a BEC of the Bloch form for D = 2 since the l/fc 2 -singularity in the Bogoliubov inequality 
(10), owing to the fc 2 -behavior of the Bogoliubov commutator (14), cannot be removed and 
still preserve the Bloch form for the BEC. For local interparticle potentials, the Bogoliubov 
commutator is given solely by the kinetic energy term of the Hamiltonian (1). This need 
not be so for nonlocal, long-range interparticle potentials {jj]. 

The fc 2 -behavior of the double commutator (14) follows from the kinetic energy term 
of the Hamiltonian since local interparticle potentials V do not contribute to the Bogoli- 
ubov commutator, viz., {\[C, V], C^\) = 0. It is interesting that the latter is not the case 
for nonlocal potentials [7]. If, for instance, the two-particle potential is a sum of a local 
and a nonlocal potential, then the former potential does not contribute to the Bogoliubov 
commutator while the latter does and if the decay of the nonlocalpotential with distance 
is sufficiently slow, then ([[C, V], &]) oc k 2 ~ e with e > as A; — )• [TJ. Accordingly, the be- 
havior of the Bogoliubov commutator in the limit k — > is dominated by the nonlocal term 
rather than the kinetic energy term as is the case in the local case. Therefore, the symmetry 
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breaking term (8) allows a BEC of the Bloch form for D = 2 provided the potential between 
the condensate atoms are given by an infinitely long-range nonlocal potential. 

Recently, a Dicke quantum phase transition was realized in an open system formed by 
a BEC coupled to an optical cavity that gives rise to a self-organized supersolid phase 
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4|. It is interesting that the phase transition is driven by infinitely long-range interactions 
between the condensed atoms. The analogy of that work to the Dicke model is based on the 
interaction Hamiltonian that gives rise to a coupling of the pump and cavity fields to the 
zero-momentum states of the atoms to the symmetric superposition of atomic states that 
carry an additional unit of photon momentum. This is quite analogous to our dynamically 
generated symmetry breaking term that allows condensation in atomic states that are integer 
multiples of a given condensate momentum. Note, however, that only one additional unit of 
photon momentum along the x and z directions are consider [4J. Higher-order momentum 
states must be included in order to describe atomic localization at t 
checkerboard pattern, that is, to represent a true BEC Bloch form 

V. SUMMARY AND CONCLUSION 



re sites of the emergent 

a. 



We use the Bogoliubov inequality to establish necessary conditions for the existence of 
two-dimensional BECs of the Bloch form and show that the interparticle potential must 
be of sufficiently long range and nonlocal. It has been shown experimentally that such a 
quantum phase transition occurs in an almost pure 87 Rb gas resulting in the onset of self- 
organization into a supersolid [4]. This behavior is equivalent to a dynamical version of the 
normal-to-superradiant quantum phase transition of the Dicke model {4]. 
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